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The problem of superluminal diffusion of relativistic particles and 
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ABSTRACT 

We discuss the problem of superluminal propagation in diffusion of ultra-high 
energy protons with energy losses taken into account. The phenomenological 



CN \ solution of this problem is found with help of generalized Jiittner propagator, 

originally proposed for the relativization of Maxwellian gas distribution. It is 
demonstrated that the generalized Jiittner propagator gives correct expressions 
in limits of diffusive and rectilinear propagation of the charged particles in mag- 
netic fields, together with the intermediate regime, in all cases without superlu- 



minal velocities. This solution, very general for the diffusion, is considered for 

two particular cases: the diffusion inside stationary objects, like e.g. galaxies, 

clusters of galaxies etc, and for the expanding universe. The comparison with 

' previously obtained solutions for propagation of UHE protons in magnetic fields 

\Q ■ 

00 



is performed. 



1 Subject headings: diffusion, magnetic fields 

O ' 
OO 
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1. Introduction 

X 

Diffusion equation describes many physical processes: thermal conductivity, the Brow- 
nian motion, the viscous motion in gases, propagation of charged particles in magnetic field 
and others. Our main interest is the diffusion of charged relativistic particles (cosmic rays). 
For collisional processes in gas, the diffusion equation can be obtained from kinetic equations 



( ILandau fc Lifshitz. v.lOl l200ll ) . The d erivation of diffusion e quation for motion of charged 



particles in magnetic fields is given in ( iBerezinsky et a/.lll990f ) 



In the extreme case, neglecting energy losses of particles and convection, assuming the 
diffusion coefficient D(t,r) independent of the time and space coordinates, the diffusion 
equation for the space density of particles n(f,t) has a simple form: 

^n(f, t) - DV 2 n(f, t) = Q(f, t), (1) 
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where Q(f,t) is the generation function. The diffusion coefficient is defined phenomenologi- 
cally through the flux density of particles as 

j(r,t) = -DVn(r,t). (2) 

Both diffusion equation ([T]) and definition of the flux density do not know about light speed 
c and in fact result in superluminal motion. A common example demonstrating it (see e.g. 



( IDunkel et all 120071 )) can be described as follows. From Eq. (T5]) one can find the streaming 



velocity u(f,t), determined from j = nu, as 

u = -DVn/n. (3) 

For the stationary (dn/dt = 0) spherically symmetric diffusion one has from Eq. (CEJ) 

and finally, using the definition of the diffusion length Id from D ~ eld for diffusion of 
relativistic particles with velocities v ~ c, one obtains from ([3]) and (J4]) 

u ~ c — , (5) 
r 

which at r < Id results in u > c. 

Another way to demonstrate the problem of superluminal propagation in diffusion is 
given by consideration of the average displacement of a particle in the diffusive regime, 
r 2 ~ Dt. One obtains for average velocity of displacement v ~ r/t ~ c(ld/r) with the same 
problem as above. From now on we shall refer to this problem as superluminal diffusion. 

From the practical point of view the described problem of superluminal diffusion, widely 
discussed in the literature, can be in many cases easily avoided. In particular, for diffusion of 
charged particles in magnetic fields, a particle deflects at length Id on average by angle 9 ~ 1 
and the movement of a particle at r < Id can be considered as quasi-rectilinear as long; as 



the e nergy spectra are concerned and as we did in our works (IBerezinsky fc Gazizov 112006 



20071 ). In these works, discussing the diffusion of ultra-high energy cosmic ray (UHECR) 
protons, we have met a more serious problem of superluminal propagation connected with 
inclusion of proton energy losses. 

The diffusion equation for protons with a source at f = r g , energy losses in the form 
dE/dt = —b(E) and with time-independent diffusion coefficient D(E), reads 

^n p (E,f,t)-D(E)V 2 n p (E,f,t) - A [b(E)n p (E, f, t)\ = Q(E,t)8 3 (f-f g ). (6) 
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The solution of this equation was found by Syrovatsky (jSyrovatsky Ill959l ) as 



1 f exp[-(r--q,)V4A(fl,£) 5 )] 

E 

where n p is the space density of relativistic protons, E is the observed energy, E g is the 
generation energy of a proton in the source, and \(E, E g ) given by 

X(E,E g )= E fdE' D{E ' ) 



E 



b(E') ' 



has a meaning of a distance squared traversed by a proton in a given direction during the 
time, when the proton energy decreases from E g to E. 

The superluminal diffusion in solution (JZJ) is immediately seen from the fact that it is the 
solution of Eq. (JBJ) only when the lower limit of integration is the observed energy i?™ 111 = E. 
One may prove it substituting ((Tj) into (jSJ). In particular, at any other lower limit E™ 11 the 
delta-function in the r.h.s. of Eq. (jHJ) is not reproduced. 

Thus, the minimum generation energy at distance r in the solution of Eq. (jHJ) equals to 
the observed energy E and it can be interpreted in the unique way: the propagation time 
t — > and thus the energy loss of a particle is absent. It implies the velocity of a particle 
v — > oo. 

The physical value of the minimal generation energy E g nin {E ) r) is given by the rectilinear 
propagation and it can be easily calculated if b(E) = —dE/dt is known. But with this _E™ in 
Eq. ([7D is not any more the solution of Eq. (jSJ). In the left panel of Fig. [T] the unphysical 
region of the solution (J7J) with superluminal velocity is shown as the hatched area. 



2. The problem of superluminal diffusive propagation of UHECR and how it 

has been eliminated 

From the formal mathematical point of view the superluminal propagation exists always 
in the solutions of the non-relativistic diffusion equations, but very often the contribution 
of unphysical regions to the solution is negligibly small (e.g. in heat conductivity). In these 
cases one may say that diffusion equation gives the correct description of the considered 
physical phenomenon. It is similar to validity of the Maxwell distribution for molecules in 
low-temperature gas (this analogy will be considered in the next section), where formally 
there are particles with velocities v — > oo. 
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Fig. 1.— The superluminal propagation problem in the Syrovatsky solution caused by 
energy losses of UHE protons for B c = 100 nG and distance to the source r = 30 Mpc. In 
the left panel is shown the region in (E g , E) plane (above E g = E line) allowed by Eq. ((7j). 
The line E g = E r g ect (E,r) corresponds to rectilinear propagation and to the physical lower 
limit in the solution (JTJ). The hatched region between E g = E and E g = E v g cct corresponds 
to superluminal velocities. For E < 1 x 10 20 eV the solutions practically do not have 
superluminal velocities. In the right panel the integrand of Eq. (jTJ) is shown as function of 
Eg for fixed E with values indicated by the numbers, given in eV. The dashed lines show the 
rectilinear physical lower limits £/™ m . At E < 1 x 10 20 eV the regions below these limits, i.e. 
with v > c are negligibly small, while at E > 1 x 10 20 eV they give considerable contribution 
toEq. ©. 
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In our works (jAloisio fc Berezinskv 1120 04. 2 0051) on diffu sion of UHE protons in static 
universe and in the works (jBerezinsky &: Gazizov 1120061 . 120071 ) on diffusion in the expanding 
universe we recognized the problem of superluminal propagation and eliminated it using the 
following recipes: 

(i) For calculation of flux from a source at distance r < ld{E) for protons with energies E 
we used the rectilinear propagation. Then the problem of superluminal diffusion disappears, 
but a technical problem of sewing the rectilinear and diffusive solutions arises, since the 
rectilinear pr opagation starts at r ~ ld( E), while diffusive regime begins at r ~ 6ld{E). 
In the paper (jAloisio fc Berezinskv 1120051) we used some s mooth interpolation between two 
regimes, while in the paper (IBerezinsky fc Gazizov 1 120071 ) we used for transition the fixed 
energy E t (r), where rectilinear and diffusive fluxes are equal. It was done intentionally, 
because the energy of transition appeared in the calculated diffuse spectrum as a feature, 
which indicates the change of the propagation regime. 

(ii) The superluminal propagation is caused also by energy losses of particles during 
diffusion, and it appears as E™ m = E in Eq. (J7J). This problem is ameliorated by presence 
of the magnetic horizon (see below): the problem of superluminal diffusion is most severe 
for the sources at largest distances beyond the magnetic horizon, while contribution of these 
sources to the observed diffuse flux is small. 

According to our calculations the picture for propagation of UHE protons in extragalac- 
tic magnetic fields looks as follows. At low energies diffusion dominates. At high energies, 
determined by energy dependence of the diffusion length ld(E) and by size of horizon rh or (E), 
the propagation changes to (quasi)rectilinear, which has no problem with superluminal prop- 
agation. This picture is valid for magnetic fields with the critical values B c (see below) 
between 0.01 — 10 nG. For the fields B c > 100 nG the superluminal signal appears at the 
highest energies in the diffuse spectrum (see Fig. 1). The fields B c < 1.6 x 10 -3 nG provide 
quasi- rectilinear propagation for all proton s at energies E > 1 x 10 17 eV, while at lower 
energies the diffusion becomes unavoidable (IBerezinsky fc Gazizov 1 120071 ). For a reasonably 
low field B c « 0.01 nG, diffusion becomes valid at E < 1 x 10 18 eV. 

Some technical details and numerical results on the problem of superluminal diffusion 
are in order. 

We have considered diffusion of UHECR in magnetized turbulent plasma with random 
magnetic field, and with the coherent magnetic field B c on the basic linear scale l c . It 
determines the critical energy E c « 1 x 10 18 (S C /1 nG)(/ c /l Mpc) eV, calculated from relation 
tl{Ec) pd l c , where tl(E) is the Larmor radius of particle in magnetic field B c . The diffusion 
length ld(E) = l c (E/E c ) 2 at E > E c for all spectra of turbulence and U{E) = l c (E/E c ) n at 
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E < E c with n = 1/3 and n = 1 for the Kolmogorov and Bohm diffusion, respectively. At 
intermediate energies the interpolation of ld{E) between these two regimes was used. For 
low magnetic field B c < 1 nG and at E > E c the diffusion length ld{E) becomes very large 
at high energy and propagation of UHE protons becomes rectilinear at all highest energies. 
The superluminal propagation thus disappears at all energies above the energy of transition 
E tr from diffusive to rectilinear regime. In calcu lation of E tr the magnetic horizon is involved 
(lAloisio fc Berezinskv 1120051 ; lAloisio et a/.ll2007l ) (for physical discussion of magnetic horizon 
see (IParizot 112004 )). Magnetic horizon is determined as a distance traversed by a particle 
during the age of a universe to: 



to 



hor 



dt D[Eg(E,t)\, 



(9) 



where E g (E, t) is the energy that a particle has at time t, if its energy at t — to is E. Putting 
dt = -dE g /b(E g ) in Eq. Q> we obtain 

-^rnax 



hor 



dE, 



D(E 9 ) 
b(E 9 ) 



(10) 



where E max = mm[E g (E, t ), E^ x ], and -E^ax is the maximum acceleration energy. 

In r^ or one may recognize the Syrovatsky variable X(E, E g ) at E g = E max . The Sy- 
rovatsky solution ([7]) knows about the horizon suppression of n(E, r) at distance r > 
a/ X(E, Eg) by exponential factor exp(— r 2 /4A). 

The energy of transition from diffusive to rectilinear propagation is determined for diffuse 
fluxes by condition ld(E) > rh or (-E), or /^(-Etr) ~ X(E tr , -E max ) with E mSuX defined as above. 



If magnetic field is very strong, B c ~ 100 nG ( lAloisio fc Berezinskv 1 120041 ). and hence 
10 20 eV, the diffusive regime continues up to very high energies and contribution of 
unphysical solutions (with v > c) appears. In Fig. [1] the integrand of Eq. ((?]) is shown as 
a function of E g for B c = 100 nG and distance to a source r = 30 Mpc. The unphysical 
regions with v > c are negligibly small at E < 3 x 10 19 eV, but is large at E = 3 x 10 20 
eV. In this case we limited our consideration by energies E < 1.2 x 10 20 eV (see Fig. 6 in 
(lAloisio & Berezinskv 112004 )). 



E c ~ 



3. Approach to solution of the superluminal propagation problem 



The most radical way to eliminate the superluminal signal in non-relativistic diffusion 
equation is given by relativization of the diffusion equation similar to relativization of the 
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Schrodinger equation in quantum mechanics. However, afte r more than 70 yea rs of efforts 
in this direction no successful solution has been found (see (IDunkel et all 120071 ) for review 
and references) . In this paper we shall follow another approach first suggested by F. Jiittner 
(jJiittner Ill91ll ) for the relativization of the Maxwellian distribution of particles. We apply 
this method to the Green function of the diffusion equation. 



In the work (IDunkel et all 120071 ) it was observed that the problem of relativization 
of the Maxwell distribution is identical to the relativization of diffusion propagator (the 
Green function). The normalized (per unit phase volume) probability density function of 
the Maxwell distribution of particles with mass m and temperature T is given by 



Pm(v) 



/ m \ 3 / 2 
\2nkTJ 



exp 



mv 
2kf 



Changing v — > x and kT/m — > 2Dt one obtains the Green function of the diffusion 
equation (CQ): 



diff 



[r,t) 



exp 



(47rDt) 3 / 2 V 4:Dt / 
where r is a distance to the source and D is time-independent diffusion coefficient. 



(12) 



Therefore, one can use for the diffusive Green function the Jiittner distribution, where 
superluminal velocities are absent. 

We start with the phenomenological definition of the propagation function (propagator). 
Let us consider first the static universe or more generically an extended stationary object 
like a galaxy or a cluster of galaxies. The propagation of the charged particles in magnetic 
fields may be described with help of propagator P(E,t,r), where E is the observed energy, 
t is the propagation time and r is the distance to a source. The definition of the propagator 
is given by 



n(E, r) 



dE n 



dtQ[E g (E,t),t] P(E,t,r) -jg(E,t) 



(13) 



where n(E, r) is the observed space density of pa rticles, E g (E,t) is energy of a particle a t 
time t, analytic expression for dE g /dE is given in (IBerezinsky. Gazizov fc Grigorieva 1120061 ). 
and Q(E,t) is the source generation function. 

The propagator P(E,t,r) can be thought of as a Green function of an unknown rela- 
tivistic equation. 

For rectilinear propagation of ultrarelativistic particles with d « c the propagator is 
given by 

^^=4^-^ (14 > 
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and Eq. (fT3j) results in 



n(E, r) 



Q[E g (E,r/c),r/c] dE g 



Ancr 2 dE ' 

which coincides with the expression obtained from conservation of number of particles. 



(15) 



For diffusive propagation the propagator can be obtained from the Syrovatsky solution 
Q as 



where 



P(r,t) 



\(E,t) 



1 



[47rA(£,t)] 3 / 2 



exp 



4X(E,t) 



rE a {t) D(E') 



(16) 



(17) 



and D(E) is the diffusion coefficient; b(E) = —dE/dt is the energy loss of a particle. The 
density of particles ndis(E,r) is given by Eq. f fl~3lh 



Both propagation functions P rec (E,t,r) and P,ns(E, t, r) are normalized by unity 



J dVP{E,t,r) = 1 



and thus they have a meaning of probability to find a particle in a unit volume at distance 
r from a source at time t after emission. It is easy to see from Eqs. f|T4l) and f|T5l) that 
P(E,t,r) has the correct dimension [P] = L~ 3 . 

The Jiittner distribution in term s of r = v and 2Dt = kT/m, as was described above, 
has been given in (IDunkel et a/.ll2006l ) as 



Pj(E,t,r) 



6{ct 



(^(§) 1-0;) 



exp 



eft 
2D 



(19) 



where 



Z(y) = 4nK 1 (y)/y (20) 

with Ki(y) being the modified Bessel function. One may observe that the superluminal 
propagation with v = r /t > c is forbidden for this propagator. 



Instead of the Jiittner function ( 1191) we introduce the generalized Jiittner function 
P g j(E,t,r), imposing to it two limiting conditions of transition to rectilinear propagator 



(TH|) and the Syrovatsky propagator f|T2l . and keeping the condition of subluminal velocities 
r < ct. For this purpose we transform the Eq. ( JT9|) as follows: 



cH 



— >■ 



cH 2 



2D 2X[E g (E,t)] 



a{t), 



(21) 



and use instead of t the new variable 



m = r/ct, (22) 

where X[E g (E,t)} is given by Eq. (fl7|) . Note that the both new quantities are dimensionless. 

The generalization imposed by Eq. ( 12 ip is motivated by the time- dependent diffusion 
coefficient D[E g (t)] and by the presence of energy losses. In this case we generalize the 
implicit quantity D x t in the Jiittner distribution f fl9|) to J D{t)dt = X(E,t) given by 
Eq. f|T7|) . As a result we have 



ID 



2+2 



C Z t 

2Dt 



2+2 



— >• 



cH 



2+2 



(ft 



a(E,t). 



(23) 



2 / D(E,t)dt 2X(E,t) 

In terms of £ and a(i?, t) the generalized Jiittner function P g j(E, t, x) and density of particles 
n(E,r) are given by: 



P 9 j(E,t,r) 



Air{ctf (1-a 2 #iK#,0] 



exp 



(24) 



n(£/, r) 



47rcr 2 



£min 



ifi[a(E,e)] 



exp 



d£ 9 



(25) 



Now we prove that Eqs. f l24l) and ( 1251) have the correct rectilinear and diffusive asymp- 
totic behavior. 

Consider first the high energy regime, where rectilinear propagation is expected. At 
E' > E c the diffusion coefficien t increases with E' as D(E ') oc (E'/E c ) 2 , while E'(t) in- 
creases exponentially with time ( lAloisio &: Berezinsky 1 120041 ) . As a result X(E,t) increases 
exponentially with time too, providing small a from Eq. (123]) and (quasi)rectilinear prop- 
agation. On the other hand from Eq. fl22|) we have £(t) ~ 1 as the second condition of 
thus (quasi)rectilinear propagation. We impose both of these conditions to Eq. ( }2"5"|) . Using 
£min = 1 — £ with e C 1 and 



a/K\{a) a 2 



valid for a -C 1, one obtains 



n(E, r) 



A 



Q 



E n 



dE g 

He 



(26) 



(27) 



with 
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where y = a/ a/1 — £ 2 , y m i n — > with a — > 0. Thus, the regime fl27|) is indeed rectilinear. 



Now consider the low-energy regime in Eq. ( f24"l) . At small E and r 3> ldis{E) the 
diffusive regime, free of the superluminal-signal problem, is valid, and we demonstrate that 
the generalized Jiittner propagator reproduces the Syrovatsky propagator, given by Eq. (J7J). 
Let us first determine the range of parameters £ and a which corresponds to the diffusion, 
described by the Syrovatsky solution. The low E corresponds to small X(E,t), and large 
enough t provides r 3> ldis(E). Therefore we have a = ct/y/X 3> 1. On the other hand 
\J X(E, t) > r provides the considerable contribution to the observed flux from a source at 
distance r and thus we have 

r 1 



< 

ct ct 



< 1. 



(29) 



With these restrictions on a and £ we compute now P g j(E,t,r) using Eq. (12^1) . At 
a 1 one has 

a/K^a) w J-a 3/2 e a . (30) 



Using exp(— a/ a/1 — £ 2 ) ~ exp[— a(l + £ 2 /2)] we obtain from ( !24j) after simple calculations 

6{ct-r) 



P(E,t,r) 



exp 



[47rA(£,t)] 3/2 

which coincides with the Syrovatsky propagator (TTdD 



4A(E,t)_ 



(31) 



However, it should be demonstrated that the space density of particles is also reproduced 
correctly. Putting P(E, t, r) from (|3TT) into ([TBI we obtain indeed the Syrovatsky solution ([7]), 
but with E™ m (r) as lower limit of integration instead of E in the Syrovatsky solution, where 
E™ m (r) corresponds to the rectilinear propagation ct = r. At very low energy £"™ m (r) — > E 
(see e.g. Fig. [1]) and both solutions coincide. 

As a matter of fact Jiittner found two versions of relativization of the Maxwell distri- 
bution. His second propagator results in the generalized propagator given by 



P 9 j{E,t,r) 



0(1 -0 I 

Aix(ctf ( 



a 5/2 K 2 [a(E,£)] 



exp 



a(E,£) 



(32) 



The difference between calculated UHE proton fluxes for both cases is less than 0.2%. In 
the following we discuss the case of propagator 



Therefore, the proposed generalized Jiittner propagator fl24|) has the correct asymptotic 
behavior (rectilinear at high energies and diffusive at low energies). It interpolates between 
these two regimes without superluminal velocities. The recent numerical simulations have 
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confirmed the Jiittner distribution of relativistic particles (ICubero et all 120071 ). which is 
generalized in our paper for the time distribution of the charged particles propagating in 
magnetic fields. 



4. Propagation of ultra-relativistic particles in magnetic fields in expanding 

universe 

Charged ultra-relativistic particles propagating in the expanding universe scatter in 
magnetic fields. Like in the previous section, where stationary flat space was considered 
(e.g. a galaxy, or static universe), in the expanding universe there are two extreme modes 
for propagation in magnetic field: the (quasi)rectilinear in the highe st energy limit and the 



diffusive in the low-energy limit. For these limits we follow the works ( iBerezinsky fc Gazizov 



20061 . 120071 ). where cosmological effects have been properly included. 



For diffusive propagation of particles from a single source located at the comoving 
coordinate x g with a generation function Q(E,t) the diffusion equation has been derived in 
the form: 

dH -b(Et) dH +3H(t)n db ( E ^ n D ^ Y*n- Q{E ' t] 5Hx x) (33) 

where the cosmological basis (x, t) is used, n(E, t, x) is the space density of the particles, 
b(E,t) = —dE/dt gives the sum of collisional energy losses, 6 int (-E), and adiabatic energy 
losses H(t)E, H(t) is the Hubble parameter, D(E,t) is a diffusion coefficient and a(t) is the 
scaling factor of the expanding universe. 

It was shown that the solution to Eq. fl33l reduces to a quadrature. Introducing the 
analog of the Syrovatsky variable 

KE,t,t>) = }dt"^^, (34) 

v 

where £ = E g (E,t,t') is a solution of the characteristic equation, the solution to Eq. fl33|) is 
given by 



n(t ,x,E) 



r z 9 


dt 


/ dz 


Jo 


dz 



exp[ -(x-x g ) 2 /A X(E,z)} 
[A7r\(E,z)f 



where z is a redshift, to is an age of the universe corresponding to z = 0, z g is the maximum 
redshift at which a source is present, dt/dz = —l/[(l+z)H(z)}, and E g (E, z) is the generation 
energy in the source. 
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The equation describing the rectilinear propagation reads ( iBerezinsky fc Gazizov 1120061 ): 
dn , , ,_, . <9n db(E.t) ce dn Q(E,t) r n.^ 

« - K«.«)flB + mt) " - 4«r " + W) « = ^ s (1 ~ *»>■ (36) 



where e is a unit vector in the direction of propagation. The solution of Eq. (1361) is given by 

n(t ,E)= g( fr g A (37) 
v °' ' Atcc x 2 s (1 + z g ) dE 1 K ' 

where x s is the comoving distance to a source. 

To use the generalized Jiittner distribution for the expanding universe we must modify 
the parameters a(t) and £(t) introduced in Section [3j 

To provide cosmologically invariant character of the parameters it is natural to substitute 
ct by the comoving length of the particle trajectory C(^) : 

at) = r — = — r dz (ss) 

Jt a{t) H J ^(1 + 2)3 + ^' 



Now we obtain a(t) from the original quantity c t/2D in the Jiittner distribution (IDunkel et al. 



20061 ). where r is propagation time, as 



c 2 r c 2 r 2 C 2 (t) 
and use a new variable £ instead of t: 

ec*) = wc(*)- (40) 

For rectilinear propagation £(t) = 1 holds as before. 

In terms of the variables a and £ defined by Eqs. ( 1391) . (l4"0l the propagator Pj(E,t,x s ) 
and particle space density n(E,x s ) have the form similar to Eqs. (j24j) and (1251) : 

P 9J (E, t, x s) =0(l-t) exp (--^=\ . (41) 



n(E x ) - _J_ f gt^j)] ^ exp f °L_A ^£ (42) 



nun 



13 



The rectilinear case corresponds to £ ~ 1 and a < 1. Integral (1421) can be evaluated 
using et/K\{ct) ~ a 2 at a < 1 and integrating over £ in a narrow range over £ « 1. One 
obtains Eq. (JHTj) as it is expected. 

The diffusive regime corresponds to a > 1 and £ < 1. Using a/Ki(a) as given by 
Eq. (|30|) one obtains for the propagator the correct diffusion expression 



exp [-x 2 /4X(E, t)] 
[4nX{E,t)f/ 2 



^,') = t::t; J] - (43) 



Thus, the generalized Jiittner propagator (!24j) has the correct asymptotic form in recti- 
linear and diffusive regimes again, providing the interpolation between them without super- 
luminal velocities. 



5. UHECR diffuse spectra in the generalized Jiittner approximation 



In this section we calculate the diffuse spectra of UHE protons propa gating in magnetic 
fields in the expanding universe and compare them with the calculations of ( iBerezinsky fc Gazizov 
20071 ). where such spectra have been calculated for diffusive propagation with transition to 
rectilinear propagation at highest energies. The proposed generalized Jiittner approxima- 
tion allows both for diffusive and rectilinear propagation at low and high energies, respec- 
tively, and provides the interpolation betwe en them based on m athematical analogy between 
Maxwellian and diffusion distribution (see ( IDunkel et all 120061 ) and Eqs. ffTTj) . ( [12"]) ). 



We start our discussion with the universal diffuse spectrum valid for homogeneous dis- 
tribution of the sources and for all modes of propagation. This spectrum can be calculated 
from the conservation of number of particles as 



n{E)dE 



to 



dt n s Q[E g (E,t)] dE, 



<:!■ 



(44) 



^min 



where n s is the source density per u nit comoving volume. According to the propagation 
theorem (lAloisio fc Berezinsky 1120041 ) this spectrum is independent of mode of particle prop- 
agation, and spectrum calculated for any specific mode of propagation must converge to the 
universal one, when distances between sources become less than any characteristic length 
involved, e.g. the diffusion length. Therefore, the spectrum (J4"2"|) integrated over space coor- 



dinates 



n(E) 



Anx 2 s dx s n s n(E, x s ) 



(45) 
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must coincide with the universal spectrum (j44[) . This fact can be easily demonstrated using 



and relations 



Atix 2 s dx s P gJ (E,t ,x s 



dt 

— and di = -— — -. 

C C 2 a(t) 



(46) 
(47) 



The calculated universal spectrum is plotted in Fig. [2j 



We calculate now the diffuse flux using the particle density from a single source given 
by Eq. (j4"2"| and assuming the lattice distribution of the sources in the coordinate space x 
with lattice parameter (the source separation) d and a power-law generation spectrum for a 
single source, 



Qs(E) 



Qoilg ~ 2 ) 
E 2 



E 
Eq 



"7s 



(48) 



where Eq is the normalizing energy, for which we will use 1 x 10 18 eV and go nas a physical 
meaning of a source luminosity in protons with energies E > Eq, L p (> Eq). The correspond- 
ing emissivity Cq = qo/d 3 , i.e. the energy production rate in particles with E > E per unit 
comoving volume, may be used to fit the observed spectrum by the calculated one. Using 
Eq. (142]) we obtain the diffuse spectrum as 



JJE) 



Qo(lg - 2) 



El 



E 



O'C [E g {E^ s )\ 



-7b 



ft 



^min 



2 s) 2 



Ma) 



cxp 



a: 



dE n 



dE 



(49) 



For the aim of comparison we us e in calculations the same ma gnetic field and its evolu- 
tion with redshift z as in the paper (IBerezinsky fc Gazizov 1 120071 ). namely we parametrize 
the evolution of magnetic configuration (l c , B c ) as 



l c (z)=l c /(l + z), B c {z) = B c (1 + z) 



2-m 



where factor (1+z) 2 describes the diminishing of the magnetic field with time due to magnetic 
flux conservation and (1 + z)~ m due to MHD amplification of the field. The critical energy 
E c (z) found from ri{E) = l c (z) is given by 



EJz) = 0.93 x 10 18 (1 + z 



1 1— TO 



Bc_ 
1 nG 



for l c = 1 Mpc. The maximum redshift used in the calculations is z B 



2. 



In Fig. [2] we compa re the generalized Jiittner s olution with the combined diffusion- 
rectilinear solution from (IBerezinsky &: Gazizov 1 120071 ). Anticipating the present paper, the 
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Fig. 2. — Comparison of the Jiittner s olution with the combined d iffusive and rectilinear 
solution (BG dotted curves) from work (IBerezinsky fc Gazizov 1 120071 ). The left panel shows 
the case B c = 0.1 nG, and the right panel B c — 1 nG, the distance between sources d = 
50 Mpc in both cases and 7 9 = 2.7. The universal spectrum is also presented for 7 9 = 2.7. 
The features seen in the BG spectra are artifacts produced by assumption about transition 
from diffusive to rectilinear propagation (see text). These features are small: note the large 
scale on the ordinate axis. 




Fig. 3.- 



The same as in Fig. [2] for B c = 0.01 nG. 
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sewing of the two solutions, diffusive and rectilinear, has been done in ( jBerezinsky &: Gazizov 



20071 ) intentionally under over-simplified assumption that for each distance to the source x g 
the transition occurs at the fixed energy E tT (x g ), where the calculated density of parti- 
cles n(E, x) for diffusive and rectilinear propagation are equal. This was done to mark 
the energy (or energy region ) of transition in the diffuse spectrum. In our earlier paper 
(jAloisio fc Berezinsky 1 120051 ) the transition was made smoothly. In Fig. [2] one can see 
tha t transition in the generalize d Jiittner solution occurs very smoothly and the feature 
in (jBerezinsky fc Gazizov 1 120071 ) solution is useful to indicate the transition. As was ex- 
plained in section [3], the various regimes of propagation are defined mostly by values of 
parameters a, in particular, a <C 1 corresponds to rectilinear propagation and a. ^> 1 to 
diffusion. For the expanding universe a(E,t) is given by Eq. (1391) . 

In Fig. H] the evolution of a(E, z) is presented as function of redshift z, along the energy 
trajectories E g = £(E,z), where E is the observed energy. The evolution is shown for 
different E and magnetic field configuration (B c ,l c ) =(0.1 nG, 1 Mpc). 

Fig. H] illustrates how protons with different observed energies E propagate in different 
regimes (diffusive and rectilinear) and how the regimes are changing. The region with a > 10 
can be considered as diffusion, a < 0.1 as rectilinear, and 0.1 < a < 10 as intermediate. 
The protons with energy E > 1 x 10 19 eV propagate only rectilinearly, E < 3 x 10 16 eV 
only diffusively, and with other energies - changing a regime of propagation. Equation for 
particle space density (J32J) takes into account automatically this changing of the propagation 
regimes. 



6. Conclusion 

Diffusion equations are intrinsically non-relativistic and superluminal velocities appear 
naturally there. The cardinal solution of this problem - the relativistic generalization of 
the diffusion equation - still expects to be found after more than 70 years of unsuccessful 
attempts. 

The phenomenological approach in its most general form consists in transition of the 
diffusive regime to the rectilinear one in all cases where superluminal propagation appears. 
We study here this problem for diffusion of relativistic particles in magnetic fields. This 
transition has a problem: there is no extreme limit in which solution of diffusion equation, e.g. 
in the form of Eq. ([7]) or Eq. fl35|) . obtains the form of rectilinear propagation. It requires an 
intermediate mode of propagation between two considered extreme regimes. This in principle 
follows from the numerical simulations. Diffusive regime is one where number of particles 
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arriving from the different directions are equal with great precision and the resulting flux is 
given by j = —DVn. In the rectilinear propagation all particles arrive from the direction 
to a source. It is clear that in intermediate regime the particles are arriving from different 
direction with noticeable asymmetry which increases as energy rises, and analytic solution 
for such regime does not exists. 

Although there is no way to obtain the rectilinear propagation as an extreme form of 
solution to diffusion equation, we do know that at high energy diffusion passes to rectilinear 
propagation. It follows from unlimited increase of D(E) with E. For example, in the case 
of diffusion in turbulent magnetized plasma with maximum scale l c , D(E) oc E 2 at E > E c 
(see section [2]). Physically it is clear that very large diffusion coefficient means rectilinear 
propagation, but this argument says nothing about possibility of the formal transition of the 
diffusive solution to the rectilinear solution. 

Our phenomenological approach for arbitrary mode of propagation is based on the 
definition of the propagator in the form of Eq. (|T3|) : 

/°° HE 
dtQ[E g (E,t),t] P(E,t,r) ~^{E,t), 

where n(E, r) is density of particles at distance r from a source with the rate of particle 
generation Q(E g ). In principle the propagator P(E, t, r) must be found as the Gree n function 



of the relativistic diffusion equation, but following the work (jDunkel et a/.ll2006l ) we found 
it from the mathematical analogy with the Maxwell distribution, given by Eqs. (lllj) and 
(02]). Then the propagator P(E, t, r) is found as relativistic Jiittner propagator, obtained for 
relativization of the Maxwell distribution, and expressed in the terms of diffusion variables 
JDunkel et artbood ). We generalized the parameters of the Jiittner distribution f|T9|) to meet 



the more general requirements of particle propagation: arbitrary dependence of diffusion 
coefficient on energy and time and including the energy losses of the particles. We proved 
that our generalized Jiittner propagator P gJ (E,t,r) given by Eqs. ( |2"lj) and (JUJ) satisfies the 
following conditions: 



• provides in all regimes velocities v < c, 

• gives in extreme cases (in particular at low and high energies) the diffusive and recti- 
linear propagators, 

• provides the normalized probability to find a particle in a unit volume at distance r 
from a source at time t after emission. Equality of this probability integrated over 
space to 1 guarantees the conservation of number of particles, 
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gives the universal spectrum for homogen eous distribution of the sou rces, as it must 
be according to the propagation theorem (lAloisio &: Berezinsky 1 12004] ) . 



In other words, following in principle (iDunkel et al we found the propagator P g j(E, t, r) 

from the general requirements listed above, which include however as the main element 
the Jiittner propagator found from the mathe matical analogy be tween the Maxwell non- 
relativistic distribution and diffusion equation ( IDunkel et a/J 120071 ). This strategy reminds 
the axiomatic approach to finding S-matrix in quantum theory, which is built on the basis 
of axioms without equations of propagation. Thus, the generalized Jiittner propagator gives 
much more than a simple interpolation between the diffusive and rectilinear regimes of prop- 
agation in magnetic fields. As one may see from the results, it successfully describes many 
features of the unified propagation and may be thought of as a reasonable approximation 
to the Green function of unknown fully relativistic equation for particle propagation, which 
includes the diffusion clS cL low-energy limit. 

We have described in this paper the calculation of UHE proton fluxes in two cases. One 
is valid for propagation in stationary objects like e.g. galaxies and clusters of galaxies. This 
solution is described in section [31 The flux from a single source at a distance r from an 
observer is given by Eq. (I25p in terms of the new variable £(£) and parameter ai(E,£). The 
diffuse flux can be calculated by summation over the sources located in vertices of the space 
lattice with the source separation d. 

The second case is given by propagation of UHE protons in the expanding universe. 
The diffuse flux is given by Eq. (jl9"|) . where variable and parameter a(E,t) are defined 
by the comoving length of a proton trajectory £(t) (see Eq. [39]) - In these calculations one 
does not need to take care of the propagation regime: it is selected automatically by values 
of £ and a in the process of integration. 
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Fig. 4. — Parameter a(E, z) as a function of redshift z for the magnetic field configuration 
(B c ,l c ) = (0.1 nG, 1 Mpc). The observed energies in EeV are shown at the evolutionary 
curves. Parameter a changes along the particle energy trajectory E g = S(E, z). The mode of 
space propagation (diffusive, intermediate, rectilinear) is changing accordingly. The regions 
of propagation correspond to diffusion (a > 10), rectilinear (a < 0.1) and intermediate 
(0.1 < a < 10). Along each trajectory the energy E g increases and a typically decreases due 
to increasing of X(E,E g ). 



